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We describe non-flat standard Friedmann cosmology of canonical scalar field with barotropic 
' fluid in form of non-linear Schrodinger-type (NLS) formulation in which all cosmological dynamical 

' quantities are expressed in term of Schrodinger quantities as similar to those in time-independent 

, quantum mechanics. We assume the expansion to be superfast, i.e. phantom expansion. We report 

■ all Schrodinger-analogous quantities to scalar field cosmology. Effective equation of state coefficient 

is analyzed and illustrated. We show that in a non-flat universe, there is no fixed WcB value for the 
j^', phantom divide. In a non-fiat universe, even w^s > —1, the expansion can be phantom. Moreover, 

in open universe, phantom expansion can happen even with WcS > 0. We also report scalar field 
exact solutions within frameworks of the Friedmann formulation and the NLS formulation in non-flat 



oo 



cr. 



X 



universe cases. 



PACS numbers; 98.80.Cq 



I. INTRODUCTION 



' Supernovae Type la data and cosmic microwave background observations show recently strong evidence of present 
^ , accelerating phase of the universe S S S IE] while nowadays inflationary paradigm in the early universe is one 
' of the corner stones in cosmology [y| . Present acceleration and inflation of the universe are both believed to result 
from effect of either dynamical scalar fleld with time-dependent equation of state coefScient Wff,{t) < —1/3 or a 
C cosmological constant with w = —I. Alternative explanation for present acceleration to dark energy is modiflcation 

of general relativity which includes braneworld models (for review, see and references therein). Among these ideas, 
■ the scalar fleld catches most attention therefore many analysis in cosmological contexts and observations have been 
I carried out Conventional formulation of canonical scalar field cosmology with barotropic perfect fiuid, can also be 
OO , expressed as non-linear Ermakov-Pinney equation as shown recently |5, fTo|. However, non-Ermakov-Pinney equation 
' for such system can also be written in form of a non-linear Schrodinger-type equation (NLS). The solutions of the 
NLS-type equation correspond to solutions of the generalized Ermakov-Pinney equation of scalar field cosmology 
[13, The NLS-type formulation was concluded and shown in case of power-law expansion in Ref. [T5| where 
. , all Schrodinger-type quantities corresponding to scalar field cosmology are worked out. NLS-type formulation also 
' provides an alternative way of solving for the scalar field exact solutions in various cases even with non-zero curvature 

til- 

Various observations allow scalar field equation of state coefficient, to be less than -1 15:]. Previous evidence 
from combined cosmic microwave background, large scale structure survey and supernovae type la without assuming 
flat universe yields = — l.OG^Q^g [16] while using supernovae data alone assuming flat universe yields = 
— 1.07 ± 0.09 [l3|. The most recent WMAP flve-year result [H, [H combined with Baryon Acoustic Oscillation 
(BAG) of large scale structure survey: SDSS and 2dFGRS[20'| and type la supernovae data from HST[31, SNLSjl 
and ESSENGE[i3 assuming dynamical w with flat universe yields —1.33 < Wd, n < —0.79 at 95% confident level [2l|. 
Also this data with additional BBN constraint of limit of expansion rate [22, l23| yields —1.29 < w^.o < —0.79 at 
95% confident level and w^.o = —1.04 ± 0.13 at 68% confident level [2l|. This suggests that the scalar fleld could be 
phantom, i.e w^j, < — 1 [24] . For a canonical scalar field, phantom behavior can be attained by negative kinetic energy 
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term of the scalar field Lagrangian density. In FLRW general relativistic cosmology, there is a Big Rip singul arity 
with a,/5, IpI — > (X) at finite time 25 1, nevertheless singularity avoidance has been attempted in various ways [26|. 



Extension to include phantom field case in NLS-type formulation was made in [l^. In NLS-type formulation one 
can presume any law of expansion a — a{t), e.g. power law a ~ or exponential expansion, a ^ exp(i/T) ITsj 
and works out all NLS-quantities keeping open possibility for the field to be phantom or non-phantom and non-zero 
spatial curvature. Analogous studies to the slow-roll, WKB and the Big Rip in NLS formulation were done in 
To attain accelerating expansion, one needs to have effective equation of state coefficient, Wcs < — 1/3 where 

WcS = , (f) 

Ptot 

Pj is density of barotropic fluid, p^ is density of the scalar field and ptot — P<j> + P-y- It ha-s been known in standard 
cosmology that for flat universe (fc — 0), if the expansion is a ~ t'' , then —I < Wcff < — f/3; if a ^ exp(t/T) , then 
Wcff = —1 and if a (ta — i)* , then w^s < —1. Here q = 2/[3(l -t- Wcff)], t, are finite characteristic times. In the 
last case, Wcs < — f corresponds to g < 0. 

In this work, we consider phantom expansion a ^ {tg. ~ tY NLS-type formulation with non-zero curvature 

k. We introduce cosmological system in Sec. [TTl then NLS-type formulation in Sec. IIIII The Schrodinger quantities 
for phantom expansion are presented in Sec. IIVI where we analyze value of WeS and show conditions of how much 
negative must be in order to keep the expansion phantom. We also illustrate parametric plots for Wcs with q and 
t. Scalar field exact solutions solved from both standard formulation and NLS-type formulation are given in Sec. |V] 
where we comment on both procedures of obtaining the solutions. Finally we conclude our work in Sec. IVII 

II. COSMOLOGICAL SYSTEM 

Barotropic fluid and scalar field fluid are major components in our scenario. The perfect barotropic fluid pressure 
and density p-y obey an equation of state, ^ — l)p-y — w^p-y while for scalar field, — w^p^. Total density 
and total pressure are ptot — P-y + P4> and ptot = P-y + Ptj>- For the barotropic fluid, w-^ is written in term of n. We set 
Wj = (n — 3)/3 so that n = 3(1 -I- w-y) = 37, hence — — f corresponds to n — — —1/3 to n = 2, = to 

n = 3, = 1/3 to n = 4, and it;.^ = 1 to n = 6. The conservation equation is hence 



P7 = -nHp^ , (2) 



with solution, 



= ■ (3) 

a" 

Therefore p.y = [(n — 3)/3](Z?/a") , where a is scale factor, the dot denotes time derivative, H — d/a is Hubble 
parameter and D > is a proportional constant. We consider scalar flcld that is minimally coupling to gravity with 
Lagrangian density, C = (l/2)e0^ — V{(j)) , of which e = 1 for non-phantom case and —1 for phantom case. Density 
and pressure of the field are given as 

P^ = + V{^) , p^ = - V{^) , (4) 



therefore 



p^ - 2F(0) 
Wd, = — — —. . 5) 

P0 e02^2F(0) 



The field obeys conservation equation 

3i/0 



0. (6) 



Considering Friedmann-Lemaitre- Robertson- Walker (FLRW) universe, the Friedmann equation and acceleration equa- 
tion are 

9 

0, Hi 

- = g-ptot(l + 3Weff) , (8) 
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where = 8ttG = 1/Mp, G is Newton's gravitational constant, Mp is reduced Planck mass and k is spatial curvature. 
Using Eqs. ©, ©, ® and (O, one can show that 



V{<f) 



H 



a 

H 2k 

y ^ 3^ 



nD 
3^ ' 



n-6\ D 
6 / ^ 



(9) 
(10) 



III. NLS-TYPE FORMULATION 



Non-linear Schrodinger-type formulation for canonical scalar field cosmology and barotro pic fluid was proposed by 
J. D'Ambroise and F. L. Williams [lll| and was also extended to include phantom field case[l2j|^. In the Schrodinger 
formulation, wave function u(x) is related to scale factor in cosmology as 



while Schrodinger total energy E and Schrodinger potential P{x) are linked to cosmology as 



(11) 



E 
P{x) 



12 

-a(t)"e0(i)2 , 



These quantities satisfy a non-linear Schrodinger-type equation: 

d2 



dx"^ 



l{x) + [E- P{x)] u{x) = -— 



with a mapping from t to a; is via 



such that hj, m 



a{t), 



x{t) — u{x) , 



±2 



Fix) 



dx . 



(12) 
(13) 

(14) 

(15) 

(16) 
(17) 



If P{x) ^ and n 7^ 0, inverse function of i)(x) exists as ^(x). Therefore x(t) = ip ^ o (p[t) and the scalar field 
potential, V o a^^{x) can be expressed as. 



12 'py-'4p(.) 



K^in? \dx ) K?n 



12^2 3fcu4/" 
^^E + — 



(18) 



IV. PHANTOM EXPANSION 

Expansion of the form a ~ (ia — i)' is called phantom when q < Q for a flat universe. Here in non-flat universe, 
q is considered to possess any value and the term phantom expansion also refers to expansion function of the form 
a ^ {ta, — as in the flat case. 



^ It is worth noting that Schrodinger-type equation in scalar field cosmology was previously considered in different procedure to study 
inflation and phantom field problems [27| . 
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A. NLS-type formulation for phantom expansion 



With the phantom expansion, a ~ (ta — ^6 use Eqs. ([TT|) and p6)) to relate Schrodinger wave function to 
standard cosmological quantity as 



U{x) = X{t) = (ta - t)-'"/' . 

Integrate the equation above so that the Schrodinger scale, x is related to cosmic time scale, t as 



x{t) 



where f3 = {qn — 2)/2 and xq is an integrating constant. Conversely, 

1 

t{x) = ta 



[P{x~xo)f^ ■ 

The Schrodinger wave function can be directly found from Eqs. (|19l) and (PTjl as 

u{x) = [/3(a; - xo)^^^''"-^'^ . 
For a (ta - 0"^: we can find e(/)(<)2 from Eq. 

2q 2k nD 



K^iU-ty K^t^-t^l 3(ta-i)«"' 

Using (|23p with phantom expansion in Eq. ([13]), therefore 



f (ia-tr-^ + ^(ta 



P(t) 

which can be expressed in term of x using Eq. (|2ip as 

2qn 1 



t) 



q{n-2) 



12 



Pix) 



(gn-2)2 {x-xoY 



kn 
~2 



{qn - 2){x - xo) 



2qin-2)/{qn-2) 



12 ■ 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



One might have a thought that all functions in phantom expansion case can be changed to those in power-law 
expansion case by interchanging (ia — t) t. However when (tg. — t) is differentiated, there is an extra minus sign. 
The Eq. I|25p slightly defers from that of the power-law expansion case because in the power-law case, the numerator 
of the second term is —2 instead of 2. The Schrodinger kinetic energy T is negative value of the first two terms of the 
Schrodinger potential. At last, the scalar field potential obtained from Eq. is 



V(t) 



g(3g-l) 2fc 



71 — 6 
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D 



(ta - ty^' 



(26) 



which can be checked by using a ^ (ta — tY in (fTO]) . Wave function of the NLS-formulation is found to be non- 
normalizable jT^ as seen Fig. in [1] for case of phantom expansion with various types of barotropic fluid. Here q is 
chosen to —6.666. In fiat universe q = —6.666 can be attained when Wcfi — —1.1. Fig. [5] shows P{x) plots for three 
cases of k with dust and radiation. In there xq = 1, therefore P{x) goes to negative infinity at x = 1. 



B. Analysis on effective equation of state coefficient 

The definition of effective equation of state coefficient, Wcff = {PtpWtf, + P-yW^) j ptot together with Eq. ^ and the 
results in Eqs. (|23p and (j26p in context of phantom expansion a ~ (ta — O'j can derive 

(-3g2 + 2g)(ta-t)2''-^-fc 
^^ff = 3,2(i^_,)2,-2+3fc • (27) 
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FIG. 1: Schrodinger wave function, u{x) when assuming phantom expansion. u{x) depends on only q, n and ta but does not 
depend on k. Here we set = 1.0 and q = —6.666. U k — 0, q = —6.666 corresponds to — —1.1. 



There is a locus, 

1/(29-2) 

t = - j , (28) 

where Weff becomes infinite along the locus. Hence for k — —1 the locus is t — tg, — g^i/(9^i) (in term of x, it is 
X = [2/{qn — 2)]g'^^"^^)/^(''"^) + a;o). Hence for fc = 0, the coefficient WeS is infinite at q = or t = ta- It seems 
from the equation above that Wcs does not depend on properties, n or amount of the barotropic fluid, D. Indeed 
WeS implicitly depends on D and n since time variable and q are related to p-y in the Friedmann equation. If fc = 0, 
it reduces to q — 2/3(1 + Wcs) and therefore the phantom condition Wcs < —1 implies q < as it is known. This 
corresponds to a condition, 

< -1- (l + w^)^ . (29) 
P<t> 

Therefore for a fluid with Wj > —1, is always less than —1 in a flat universe. In order to have the expansion 
a ^ (^a — t)"^ in A; = universe, we must have WcS < — 1, i.e. in phantom region. We can rewrite in term of WcS as 

[^(ia-t)"' + ^(ia-t)"'1Wcff- ^£'(ta-i)"«" , , 

^(ia - + f (ta - i)-"' - ^(ia - t)-^" 

Eq. (j30p . when D — and k — 0, yields = WcS- Albeit we set only D = 0, it gives the same result since is 
independent of geometrical background. However, since the expansion law is fixed, is tied up with D implicitly 
via Eq. ([1]). Note that has value in the range (— oo, —1] and [1, oo) so that the phantom crossing can not happen 
when the scalar field is dominant. However, presence of the dust barotropic fluid in the system gives a multiplication 
factor that is less than 1 to the equation of state, i.e. 

u>es = ( — — ) '^'P ■ (31) 

We can see that the phantom crossing from WcS > — 1 to Wcs < — 1 can happen in this situation. Fig. [3] presents 
parametric plots of the W(.f[-,q,t diagram for various k values. From the figure, we see the locus in Eq. (|28p where 
Weft blows up. In the parametric plots, the value of Wcff at any instance can be obtained if we know the value of q. 
We need to know q from observation in order to know the realistic value of Wcff or the other way around. Fig. [4] 
plotted from Eq. (I27p setting ta = 1 and t = 0.7 shows that ii k = ±1, q could be negative, i.e. phantom accelerated 
expansion, even when Wcff > — 1. Regardless of ta and t, 

lim Weff(g) = -1 and lim Weff{q) ^ -7: , (32) 

9 — ^ — 00 q — >+oo o 
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Scalar field with dust (n = 3) Scalar field with radiation (n = 4) 
P(x) Fix) 



k = 



Fix) 



P(x) 



P(x) 



P(x) 



FIG. 2: Schrodinger potential in phantom expansion case for dust and radiation fluids with fc = 0, ±1. Numerical parameters 
are as in the u{x) plots (Fig. [TJ. xo is set to 1. For non-zero k, there is only one real branch of P{x). 

for phantom expansion. In particular, for k — —1^ WcH > could give g < and Wos is infinite when In q/ In (ta ~ t) + 
q=l (see Eq. 

V. SCALAR FIELD EXACT SOLUTION 
A. Bound value of <^(f) from effective equation of state for A: = case 
In flat universe, the phantom expansion occurs when Wcfr < — 1. Using Eqs. ([5]) in Eq. ([1]), we get a bound 



£0 < -gPT 



(33) 



Assuming a{t) — {ta — t)'^ and phantom scalar field, i.e. e — —1 with using Eq. ([3]), the solution is found to be in the 
region. 



, , , 1 Dn , , ,5 



(34) 



where (3 = {qn — 2)/2 . 



FIG. 3: Parametric plots of WcS for the expansion a ~ (fa — ty in closed, flat and open universe. Here ta. is set to 1. 

B. Solution solved from Priedmann equation 

1. Scalar field potential in flat and scalar field dominated case 

A simplest case for analysis is when considering flat universe (k — 0) with negUgible amount of barotropic fluid 
{D — 0). The Eq. ([23]) is hence simply integrated out. The solution is 

(/.(t) = ±i^/^ln(ta-0 + 00 (35) 
K V e 

Insert this result into Eq. (pS)) . we obtain the scalar field potential, 

y(^) = ^^^^^ exp S^±K^m) - ^o]} . (36) 
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2 4 
_i I I I i_ 



1 

'W^[f =-1/3 



k = 





FIG. 4: w^ff for the expansion a ^ (ta — i)' in closed, fiat and open universe. Here ta is set to 1 and t is 0.7. 



The solutions above are real only when q and e have the same sign, i.e. when e = l,q > and e = —1, g < 0. This 
looks similar to potential that gives power-law expansion as well-known [2^. It is not surprised since in our case 
{q < 0) it has been known that phantom field, when rolling up the hill of slope- varying exponential potential (varying 
q), results in phantom expansion a ^ {ta — t)'^ psl ]. 
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2. Solution for k = 0, D ^ case 



For the case fc = with D ^ 0, the solution of Eq. (|23p is 



qn — 2 \ eti'' 



hi 



(l + - inDKy6q){t, ~ t)-'?"+2y 



6q 




(37) 

which is infinite when q = 2/ri. The last logarithmic term in the bracket is an integrating constant. Logarithmic 
function is valid only when q < 0. 



3. Solution for k ^ 0, D — case 



For the reverse case, fc 7^ 0, D = 0, the solution is 

g — 1 V 



In 



(ia - + 1 



(38) 

which becomes infinite when q = I. The values q and e must have the same sign for it to be real- value function. The 
case fc ^ with D ^ can not be found analytically except when setting n = 2 {w~f — —1/3) which is not natural 
fluid. 



C. Solution solved with NLS-type formulation 



One can obtain exact solution of Eq. ((23|) indirectly via NLS-type formulation. Consider Eq. ([25|) , we notice that 
setting D — does not make sense in NLS-formulation since even D vanishes, n (barotropic fluid parameter) still 
appears in other terms. Therefore we can only consider non-zero D case. Assuming fc = with D ^ and using Eq. 
([551) in Eq. the solution is 



eK'^{qn — 2)^ 



-^ 1 



K?'Dn{qn — 2)^ 
2Aq 



(x - x^f 



In 



1 + yjl - [K'^Dn{qn - 2)2/24g] (x - xqY Aqn 



(x - xo) 



e{qn - 2) 



(39) 



Transforming to t variable using Eq. (|20p . 



m = ± 



1 / 2q 
qn — 2\ en'^ 



In 



(ta - t) 



-qn+2 



(1 + ^1- inDKy6q)it,~t)-in+2y 



2WI 



nDK2(ta - i)-«"+2 /qn~2 



6q 



In 



\ 2qn 



(40) 



The only difference from the solution ([37]) obtained from standard method is the logarithmic integrating constant 
term in the bracket. In case of fc 7^ with D 0, the integral (11711 can not be integrated analytically even when 
assuming n value except for n — 2 which is integrable. However n = 2 is not natural fluid. This is similar to using 
standard method in Sec. IVB 31 
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VI. CONCLUSIONS 

We consider a system of FLRW cosmology of scalar field and barotropic fluid assuming phantom acceleration. We 
have worked out cosmological quantities in the NLS-formulation of the system for flat and non-flat curvature. The 
Schrodinger wave functions are illustrated in Fig. [T] for various types of barotropic fluid. These wave functions are 
non-normalizable. We show Schrodinger potential plots for dust and radiation cases in closed, flat and open universe. 
The procedure considered here is reverse to a problem solving in quantum mechanics in which the Schrodinger 
potential must be known before solving for wave function. In NLS formulation, the Schrodinger equation is non-linear 
(reducible to linear in some cases) and the wave function is expressed first by the expansion function, a{t). Afterward 
the Schrodinger potential is worked out based on expansion function assumed. Moreover, the NLS total energy E is 
negative (see Eq. ([12]) )• We also perform analysis on effective equation of state. We expresses WeS in term of q and 
k. In a non-flat universe, there is no fixed WcS value for a phantom divide. We show this by analyzing Eq. p7p and 
by presenting illustrations in Figs. ([3]) and In these plots, even w^s > — 1, the expansion can still be phantom, 
i.e. q can be negative. Especially, in fc = — 1 case, positive WcS could also give g < 0. The value of Wcff approaches -1 
when q — oo and —1/3 when q +oo. In open universe, w^s blows up when Ing/ In (tg^ — t) + q = 1. 

The last part of this work is to solve for scalar field exact solution for phantom expansion. Within framework 
of the standard Friedmann formulation, we obtained exact solution in simplest case where scalar field is dominated 
in flat universe. Apart from that we also obtained exact solutions in the cases of non-flat universe with scalar field 
domination and flat universe with mixture of barotropic fluid and scalar field. Afterward, we use NLS formulation, in 
which the wave function is equivalent to the scalar field exact solution, to solve for the exact solutions. We can apply 
the NLS method to solve for the solution only when the barotropic fluid density is non-negligible. Setting D — 
in NLS framework is not sensible because even if D term vanishes, the barotropic fluid parameter n still appears in 
other terms of the wave function. This is a disadvantage point of the NLS formulation. 

Transforming standard Friedmann formulation to NLS formulation renders a few effects to the integration. In 
standard form (Eq. n appears in only D-term and all terms are i-dependent. In NLS-form (Eq. (P5|) when 

inserted in Eq. ^Tf} ). D-teim becomes a constant (E), hence the number of a;(or equivalently t)-dependent terms is 
reduced by one. This is a good aspect of the NLS. In both Friedmann-form and NLS-form, the solutions when fc ^ 
and D 7^ are difficult or might be impossible to integrate unless assuming values of q and n. Therefore reduction 
number of x-dependent term helps simplifying the integration. 
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